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f— *^ , I discuss a model for quantized gravitation based on the simplicial lattice dis- 

^v^j i cretization. It has been studied in detail using a comprehensive finite size scaling 

analysis combined with renormalization group methods. The results are consistent 
G ' with a value for the universal critical exponent for gravitation v = 1/3, and sug- 

gest a simple relationship between Newton's constant, the gravitational correlation 
length and the observable average space-time curvature. Some perhaps testable 
(^-^ , phenomenological implications are discussed, such as the scale dependence of New- 

^s^j ■ ton's constant and properties of quantum curvature fluctuations. 

^ . 1 Introduction 

in 

' One of the outstanding problems in physics is a determination of the quantum- 

mechanical properties of Einstein's theory of Gravitation. Approaches based on 
linearized perturbation methods have had limited success so far, as the underlying 
theory is known not to be perturbatively renormalizable. Since gravitational fields 
are themselves the source for gravitation already at the classical level, perturbative 
results are perhaps of doubtful validity for sufficiently strong effective couplings, 
'"pi ' This is especially true in the quantum domain, where large fluctuations in the 

JL gravitational field appear at short distances. In general nonperturbative effects 

can give rise to novel behavior in quantum field theory, and in particular to the 
r"| , emergence of non-trivial fixed points of the renormalization group, a phase transi- 

tion in statistical mechanics language. Then the universal low energy behavior of 
field theories is almost completely determined by the fixed point structure of the 



o 



^ ' renormalization group trajectories 



The situation described above clearly bears some resemblance to the theory of 
strong interactions, Quantum Chromodynamics, where non- linear effects are known 
here to play an important role, and end up restricting the validity of perturbative 
calculations to the high energy, short distance regime □. For low energy prop- 
erties the lattice formulation, combined with renormalization group methods and 
computer simulations, has provided convincing evidence for quark confinement and 
chiral symmetry breaking, two phenomena which are invisible to any order in the 
usual perturbative expansion. 

2 The Lattice Theory 

A discrete lattice formulation can be applied to the problem of quantizing gravita- 
tion; instead of continuous metric fields, one deals with new gravitational degrees 
of freedom, the edges, which live only on discrete space-time points and interact 
locally with each other. In Regge's simplicial formulation of quantum gravity a 
one then approximates the functional integration over continuous metrics by a dis- 
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cretized sum over piecewise linear simplicial geometries uU. In such a model the 
role of the continuum metric is played by the edge lengths of the simplices, while 
local curvature is described by a set of deficit angles, determined via known for- 
mulae as functions of the given edge^lengths. The simplicial lattice formulation 
of gravity is locally gauge invariant and can be shown to contain perturbative 
gravitons in the lattice weak field expansion 0, making it an attractive and faithful 
lattice regularization of the continuum theory. In the end the original continuum 
theory of gravity is to be recovered as the space-time volume is made large and the 
fundamental lattice spacing of the discrete theory is sent to zero, without having to 
rely, at least in principle, on any further approximation to the original continuum 
theory. 

Quantum fluctuations in the underlying geometry are represented in the discrete 
theory by fluctuations in the edge lengths, which can be modeled by a well-defined, 
and numerically exact, stochastic process. In analogy with other field theory models 
studied by computer, calculations are usually performed in the Euclidean imaginary 
time framework, which is the only formulation amenable to a controlled numerical 
study, at least for the immediate foreseeable future. The Monte-Carlo method, 
based on the concept of importance sampling, is well suited for evaluating in a 
numerically exact way the discrete path integral for gravity and for computing the 
required averages and correlation functions. By a careful and systematic analysis of 
the lattice results, the critical exponents can be extracted, and the scaling properties 
of invariant correlation functions determined from first principles c3. 

The starting point for a non-perturbative study of quantum gravity is a suitable 
definition of the discrete Feynman path integral. In the simplicial lattice approach 
one starts from the discretized Euclidean path integral for pure gravity, with the 
squared edge lengths taken as fundamental variables, 



z L = n (^( s )) ct n di i®i i u ex p \ - e ( a Vh - k ^ + a s i A i/ v h) ■ 

J° s ij { h J 

(1) 

The above expression represents a lattice discretization of the continuum Euclidean 
path integral for pure quantum gravity 



(2) 

(with k 1 = 8ttG, and G Newton's constant) and reduces to it for smooth enough 
field configurations. The discrete gravitational measure in Zl should be considered 
as the lattice analog of the DeWitt continuum functional measure. The SA term 
in the lattice action is the well-known Regge term 0, and reduces to the Einstein- 
Hilbert action in the lattice continuum limit 0. A cosmological constant term 
is needed for convergence of the path integral for large edge lengths, while the 
curvature squared term allows one to control the fluctuations in the curvature □. 
In general in the Regge theory the correspondence between lattice and contin- 
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uum operators is 



hinges hZ)x 



^gR{x) -> 2 ^ 



hinges fOa; 



(3) 



hinges /Ox 



where (5^ is the deficit angle at the hinge (triangle) /i, is the area of hinge labelled 
by h, and is the four-volume associated with the same hinge h. 

In practice and for obvious phenomenological reasons one is only interested in 
the limit when the higher derivative contributions are small compared to the rest, 
a — > 0. In this limit the theory depends, in the absence of matter and after a 
suitable rescaling of the metric, only on one bare parameter, the dimensionless cou- 
pling k 2 /A. Without loss of generality, one can therefore set the bare cosmological 
constant A = 1. In the following discussion only the case a = will be considered. 

3 Phases of Lattice Gravity 

Studies oiLsmall lattices suggest a rich scenario for the ground state of quantum 
gravity uau. The present evidence indicates that simplicial quantum gravity in 
four dimensions exhibits a phase transition (in the bare coupling G) between two 
phases: a strong coupling phase, in which the geometry is smooth at large scales 
and quantum fluctuations in the gravitational field eventually average out and are 
bounded; and a weak coupling phase, in which the geometry is degenerate and 
space-time collapses into a lower-dimensional manifold, bearing some physical re- 
semblance to a spiky branched polymer. Only the smooth, small negative curvature 
(anti-DeSitter like) phase appears to be physically acceptable. 

Phrased in a different language, the two phases of quantized gravity can loosely 
be described as having in one phase (with bare coupling G < G Cl the branched 



with a vanishingly small negative average curvature in the vicinity of the critical 
point at G c . 

The existence of a phase transition at finite coupling G, usually associated 
in quantum field theory with the appearance of an ultraviolet fixed point of the 
renormalization group, implies in principle non-trivial, calculable non-perturbative 
scaling properties for correlations and effective coupling constants, and in particular 
here for Newton's gravitational constant. Since only the smooth phase with G > G c 
has acceptable physical properties, one would conclude on the basis of fairly general 
renormalization group arguments that at least in this lattice model the gravitational 
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(4) 



(5) 



coupling can only increase with distance. In addition the rise of the gravitational 
coupling in the infrared region rules out the applicability of perturbation theory 
to the low energy domain, to the same extent that such an approach is deemed 
to be inapplicable to study the low-energy properties of asymptotically free gauge 
theories. 

It is a remarkable property of quantum field theories that a wide variety of phys- 
ical properties can be determined from a relatively small set of universal quantities. 
Namely the universal leading critical exponents, computed in the vicinity of some 
fixed point (or fixed line) of the renormalization group equations. In the lattice 
theory the presence of a fixed point or phase transition is often inferred from the 
appearance of non-analytic terms in invariant local averages, such as for example 
the average curvature defined as 

2 <JJ^gR{xy> _ Av-l 
K 1 > <Jd*x^g> = n{k) k~k c ^ - fc ) • < 6 > 

From such averages one can infer the value for v 1 the correlation length exponent, 

m ~ A i {k c -k)- v . (?) 

k— »fc c 

In terms of the physical correlation length £ one has 

kZ ^ e /v ~ A ■ (8) 

Correct dimensions here can be restored by supplying appropriate powers of the 
ultraviolet cutoff, the Planck length lp = \[G. The fundamental critical exponent 
v is related to the derivative of the beta function for G 

BC 

/3(G) = (9) 

olog^i 

in the vicinity of the ultraviolet fixed point at G c , 

p'(G c ) = -l/v . (10) 
Integrating close to the non-trivial fixed point one obtains for G > G c 

= «"' = Aex »(-/ G sl))c~ c « A '^- i '- w ' <CJ ■ <") 



m 



where m is an integration constant identified with the inverse correlation length, 
and A — lp 1 the ultraviolet cutoff. 

This physical correlation length £ determines the long-distance decay of the 
connected invariant two-point correlations at fixed geodesic distance d. For example 
for the curvature correlation one has, for distances much larger compared to the 
correlation length £, 

<VgR{x)y/gR{y)5(\x-y\-d)> c ~ dr" e~ d ^ , (12) 
while for short distances one expects a power law decay 

<V9R(x)y/gR{y)5{\x-y\-<t)> Cd ~ i d2( l_ 1/v) ■ (13) 
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Method 


kc 


V 


TZ vs. k 


0.0630(11) 


0.330(6) 


TZ a vs. k 


0.0639(10) 




Xn vs. k 


0.0636(30) 


0.317(38) 


3/2 , 
Xii vs. k 


0.0641(17) 




X n/((l 2 )n) vs. k 


0.0635(11) 


0.339(9) 


Xn vs. TZ 




0.328(6) 


1Z FS scaling 




0.333(2) 


Xk FS scaling 




0.318(10) 



Table I: Summary table for the critical point k c and the critical exponent v , as obtained 
from lattices with up to 16 4 sites. The last three entries assume a critical point at 

k c = 0.0636. 



One further result is that according to the renormalization group the scale depen- 
dence of the effective Newton constant is given by 

G(r) = G(0) [l + c(r/0 1/u + O((r/0 2/ ^)] , (14) 

with c a calculable numerical constant. Here the scale £ _1 plays a role very similar 
to the scaling violation parameter Ajj^ of QCD. It seems natural, although para- 
doxical at first, to associate £ with some macroscopic cosmological length scale. 

It should be clear from this brief discussion that the critical exponents pro- 
vide a wealth of useful information about the continuum theory. In reality, the 
complexity of the lattice interactions and the practical need to sample many sta- 
tistically independent field configurations contributing to the path integral leads to 
e requirement of powerful computational resources. The results discussed here 
3 were obtained using a dedicated custom-built 20-GFlop 64-processor parallel 
supercomputer, described in detail elsewhere 113. 

Table I 1121 summarizes the results obtained for the critical point k c = l/8nG c 
and the critical exponent v. From the numerical calculations one finds 

k c = 0.0636(11) v = 0.335(9) , (15) 

which clearly suggests v = 1/3 for pure quantum gravity. 



4 Critical Exponents and Phenomenology 

This section contains a brief discussion of some of the consequences that follow from 
the results presented in the previous section, and in particular the result v = 1/3. 

First let us notice that the value v = 1/3 does not correspond to any known 
field theory or statistical mechanics model in four dimensions. Indeed for all scalar 
field theories (spin s = 0) in four dimensions it is known that v = 1/2, while for the 
compact Abelian U(l) gauge theory (s — 1) one has v — 2/5. The value v = 1/3 
for s — 2 is then consistent with the simple interpolation v^ 1 = 2 + s/2. 

One distinctive feature of the results is the appearance of a gravitational cor- 
relation length £. Naively one would expect, simply on the basis of dimensional 
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arguments, that the curvature scale should get determined by this correlation length 

n rZo i/e , (i6) 

but one cannot in general exclude the appearance of some non-trivial exponent. In 
the previous section arguments have been given in support of the value v = 1/3 for 
pure gravity. From the equation relating the average curvature to the gravitational 
correlation length one has 

^(0 ~ ~~o — 7TTT • ( 17 ) 

Here the correct dimension for the average curvature 1Z has been restored by sup- 
plying appropriate powers of the ultraviolet cutoff, the Planck length lp = \/~G. 
For v — 1/3 in four dimensions one then obtains the remarkably simple result 

An equivalent form can be given in terms of the curvature scale Hq, defined through 
R = —12Hq which has dimensions of a mass squared. Then close to the critical 
point 

Hi = Cniipm , (19) 

where /ip = 1/y/G is the Planck mass, m — l/£ is the inverse gravitational corre- 
lation length, and Ch ~ 4.9 a numerical constant of order one (the value for C is 
extracted from the known numerical values for TZ and rn close to the critical point 
at fc c ). It is amusing to note that this result is reminiscent of the pcac relation in 
pion physics. 

One can raise the legitimate concern of how these results are changed by quan- 
tum fluctuations of matter fields. In the presence of matter fields coupled to gravity 
one expects the value for v to change due to vacuum polarization loops containing 
these fields. But based on general arguments, one would expect fields whose masses 
are significantly above m = l/£ to give negligible contributions to vacuum polar- 
ization loops, and thus leave the universal critical exponents which characterize the 
large distance behavior unaffected. 

It seems natural to identify Hq with either some (negative) average spatial 
curvature, or possibly with the Hubble constant determining the macroscopic ex- 
pansion rate of the present universe u. In the Friedmann-Robertson- Walker model 
of standard cosmology on has for the Ricci scalar 



Rmcci = -6 - + -== + -} , (20) 




where R(t) is the FRW scale factor, and k = 0, ±1 for spatially flat, open or closed 
universes respectively. Today the Hubble constant is given by Hq = (R/R)1 but it 
is eventually expected to show some slow variation in time. Its characteristic length 
scale today cH^ 1 as 10 28 cm is comparable to the extent of the visible universe. 
Under such circumstances one would expect the gravitational correlation £ to be 
significantly larger than cH^ 1 . 
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A potential problem arises though in trying to establish a relationship between 



quantities which are truly constants (such as the ones appearing in Eq. (19)), and 
H which most likely depends on time (the only exception being the steady state 
cosmological models, where H is truly a constant of nature; these models are not fa- 
vored by present observations, including detailed features of the cosmic background 
radiation) . In any case it is clear that some of these considerations are in fact quite 
general, to the extent that they rely on general principles of the renormalization 
group and are not tied to any particular value of u, although the favored value 
v = 1/3 clearly has some aesthetic appeal. 

One further observation can be made regarding the running of G in Eq. (fTi]). 
Assuming the existence of an ultraviolet fixed point, the effective gravitational 
coupling is given for "short distances" r <C £ by 

G(r) = G(0)[l + c(r/£) 3 + O((r/0 6 )] , (21) 

with c a calculable numerical constant of order one. The appearance of £ in this 
equation, which is a very large quantity by Eq. (|l9|), suggests that the leading 
scale-dependent correction which gradually increases the strength of the effective 
gravitational interaction as one goes to larger and larger length scales, should be 
exceedingly small. It also suggests that the deviations from classical general rela- 
tivistic behavior for most physical quantities is in the end practically negligible. 

It is only for distances comparable to or larger than £ that the gravitational 
potential starts to weaken and fall off exponentially, with a range given by the 
gravitational correlation length £, 

V(r) ~ -G(r) ^ ^ . (22) 

r » £ r 

In many ways these results appear qualitatively consistent with the expected be- 
havior of the tree-level graviton propagator in anti-de Sitter space. In the real 
world the range £ must be of course very large; from the fact that super-clusters of 
galaxies apparently do form, one can set an observational lower limit £ > 10 25 cm. 

It is unclear to what extent gravitational correlations can be measured directly. 
From the definition of the curvature correlation function in Eq. ( |l3| ) one has for 
"short distances" r -C £ and for the specific value v — 1/3 the remarkably simple 
result 

<JgR{x) ^R(y)5(\x-y\-d)> c ± , (23) 

with A a calculable numerical constant of order one. One should contrast this 
behavior with the semiclassical result attained close to two dimensions (and which 
incidentally coincides with the lowest order weak field expansion result), which gives 
instead for the power the value 2(rf — 1 jv) ~ 2(d — (d — 2)) ~ 4, as expected on the 
basis of naive dimensional arguments (R ~ d 2 h). 

Next one can look at local curvature fluctuation. If one considers the curvature 
R averaged over a spherical volume V r = 47rr 3 /3, 

= it / y/g&t) R&t) (24) 
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one can compute the corresponding variance in the curvature 

[S(V9 R)f = J <?x cPy < R(x) ^ R{y) > c = |A . (25 ) 

As a result the r.m.s. fluctuation of ^fgR averaged over a spherical region of size r 
is given by 

S(V~9R) = ^ \ , (26) 
while the Fourier transform power spectrum at small k is given by 

^ = l^%| 2 = ^i- 07) 

One can use Einstein's equations to relate the local curvature to the (primordial) 
mass density. From the Einstein field equations 

Rvu, - y^R = 8ttT^ (28) 

for a perfect fluid 

T\iv = V9ixv + (P + p)u^u u (29) 
one obtains for the Ricci scalar, in the limit of negligible pressure, 

R(x) w 8ttG p{x) . (30) 
As a result one expects for the density fluctuations a power law decay of the form 

< P(x) P (y) > c 7 1 2 . (31) 

\x-y\ 2 

Similar density correlations have been estimated from observational data by ana- 
lyzing known galaxy number density distributions, giving a value for the exponent 
of about 1.77 ± 0.04 for distances in the Wkpc to lOMpc range. 

5 Outlook 

Numerical simulation methods combined with modern renormalization group ar- 
guments can provide detailed information on non-perturbativc aspects of a lattice 
model of quantum gravity. One then finds that the lattice model has two phases, 
only one of which is physically acceptable. In spite of the fact that the Euclidean 
theory becomes unstable as one approaches the critical point at k c , it is still possi- 
ble to determine by a straightforward analytic continuation the physical properties 
of the model in the vicinity of the true fixed point, defined as the point where 
a non-analiticity develops in the strong coupling branch of Z^ik). There scaling 
implies that the physical correlation £ diverges. 

If this prescription is followed, an estimate for the non-perturbative Callan- 
Symanzik beta function in the vicinity of the fixed point can be obtained, to leading 
order in the deviation of the bare coupling from its critical value. The resulting 
scale evolution for the gravitational constant is then quantitatively quite small, if 
one assumes that the scaling violation parameter is related to an average curvature 
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and its characteristic scale H . Its infrared growth, consistent with the general idea 
that gravitational vacuum polarization effects cannot exert any screening, suggests 
that low energy properties of quantum gravity are inaccessible by weak coupling 
perturbation theory: low energy quantum gravity is a strongly coupled theory. As 
pointed out in the discussion there are a number of attractive features to the pure 
gravity result v = 1/3, including a simple form for the curvature correlations at 
short distances. 

It seems legitimate to ask the question whether the present lattice model for 
quantum gravity provides any insight into the problem of the cosmological constant. 
The answer is both yes and no. To the extent that a naive prediction of quantum 
gravity is that the curvature scale should be of the same order of the Planck length, 
1Z ~ 1/G, the answer is definitely yes. Indeed it can be regarded as a non-trivial 
result of the lattice models for gravity that a region in coupling constant space can 
be found where space-time is stiff and the curvature can be made much smaller than 
1/G. In fact the evidence indicates that the average curvature 1Z vanishes at the 
critical point k c . And this is achieved with a bare cosmological constant A which 
is of order one in units of the cutoff. Phrased differently, the dimensionless ratio 
between the renormalized and the bare cosmological constant becomes arbitrarily 
small towards the critical point. 

At the same time the effective long distance cosmological constant is clearly non- 
vanishing and of order l/£ as a consequence of dimensional transmutation. The 
value zero is only obtained when £ is exactly infinite, which happens only at the 
critical point k c . Thus to make the effective cosmological constant small requires a 
fine tuning, in the sense that the bare coupling k c — k has to be small. But since the 
correlation length determines the corrections to the Newtonian potential (and in 
particular its eventual decrease for large enough distances) , it would seem unnatural 
to have a short correlation length £: in such a world there would be no long- 
range gravitational forces, and separate space-time domains would have decoupled 
fluctuations. From this perspective, long range forces and a small cosmological 
constant go hand in hand. Quantum fluctuation effects show that hyperbolic space- 
times with small curvature radii cannot sustain long-range gravitational forces, at 
least in this lattice model. 
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